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We study the stabilization of localized structures by discreteness in one-dimensional lattices of 
diffusively coupled nonlinear sites. We find that in an external driving field these structures may lose 
their stability by either relaxing to a homogeneous state or nucleating a pair of oppositely moving 
fronts. The corresponding bifurcation diagram demonstrates a cusp singularity. The obtained 
analytic results are in good quantitative agreement with numerical simulations. 
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The problem of dynamics in discrete nonlinear lattices 
arises in diverse physical, biological, and engineering sys- 
tems. Among the examples are interaction of charge- and 
spin-density waves with impurities in correlated electron 
materials arrays of Josephson junctions Q, calcium 
release waves in living cells || , systems of coupled nonlin- 
ear oscillators etc. Of special interest are lattices with 
continuous coupling between the nonlinear sites, because 
of their rich dynamical behavior. In Ref. |B| such systems 
have been shown to carry propagating burst waves. The 
role of localized structures as nucleation embryos in driven 
nonlinear systems is of wide importance ||, but the in- 
fluence of lattice discreteness on nucleation and transport 
has received much less attention to date. We will show 
in this Letter that lattice discreteness indeed results in 
qualitatively new phenomena. 

Below we study the stabilization of localized structures 
by discreteness in a system consisting of a lattice of diffu- 
sively coupled nonlinear sites. We investigate the proper- 
ties of such structures and different mechanisms of their 
instabilities. We find a hysteresis effect in the nucleation 
of propagating fronts from the localized structures in an 
external field. We obtain our results analytically and 
confirm them by numerical simulations of the full sys- 
tem dynamics. Our system has relaxational dynamics 
given by du/dt = —8£/5u , where Z is the system energy 
functional, u is the order parameter, and 5/Su is a vari- 
ational derivative. We consider the energy functional of 
the following form 

m = § /(£) 2 ^ + «E^n- (i) 

i 

Here D is the diffusive coupling, T is the discrete non- 
linear potential, a is the potential amplitude, and the 
sum is over all lattice sites i . We study the behavior of 
system ([!]) with bistable dynamics. The details of the 
potential shape are not important so long as it has at 
least two minima, u = u_ and u — u + separated by a 
barrier with a maximum at u = uq . The simplest exam- 
ple of such a potential T is a fourth degree polynomial 
(4> 4 ). Here we study the sine-Gordon potential, because 



of its applicability to modeling the dynamics of charge 
density waves [|J . Although this potential has an infinite 
number of wells, for the localization problem studied in 
this Letter, only two neighboring minima are relevant. 
The potential has a form 

J-[ui] = — cosuj — Eui, (2) 

where E is the applied external field. The field is required 
to make the depths of the potential wells different and 
thus allow for the propagation of fronts from less stable to 
more stable states. For potential (||), these stable states 
are it_ = arcsini? , u + = u_ + 2tt , uq = tt — u_ . The 
system dynamics are then governed by 

u t = (3u xx + )^5(x-i)(-sinu + E) , (3) 

i 

with the space and time rescaled as x — ► xjdid is the dis- 
tance between sites) and t — ► ta , respectively. The renor- 
malized dimensionless diffusive coupling is (3 = D/ad 2 . 

When (3 1 , the system is in the continuous limit, 
where it is described by the overdamped sine-Gordon 
equation: u t = (3u xx — sinu + E . If one starts with 
the bulk of the system in the state u = U- , below the 
barrier maximum uq , and a finite size nucleus of the state 
u = u + , above the maximum, then, for E > , there ex- 
ists a critical size of the nucleus. If the nucleus exceeds 
this critical size, it breaks into a pair of oppositely mov- 
ing fronts; otherwise it relaxes back to the bulk state. 
We show that introducing discreteness into the system 
can stabilize the critical nucleus as localized structure. 

In the discrete regime, j3 < 1 , the system demon- 
strates substantially richer behavior. Well-separated 
fronts undergo a ■pinning- depinning transition from sta- 
tionary kinks to propagating burst waves, which are peri- 
odic in a traveling- wave reference frame j5| . A stabilized 
nucleus represents a bound pair of kink and antikink. 
There exists a hierarchy of nuclei with differing numbers 
of sites at the state u = u + , above the maximum. The 
distance between the kink and the antikink increases with 
the number of these sites, and the binding energy rapidly 
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decays. Since the nucleation is caused by local fluctua- 
tions, the physically most important nuclei are those with 
small numbers of sites. As the parameters vary, a nucleus 
may lose its stability. Fig. |l| shows two alternative sce- 
naria of one-site nucleus destabilization: (a) breaking to 
a pair of burst waves and (b) relaxation to the bulk state. 



equations for iij's, the values of u at the sites x^s, 
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FIG. 1. Successive profiles of the numerical solution u(x) 
of Eq. (|^) , starting from a one-site nucleus, (a) Nucleation of 
two oppositely moving burst waves; E — 0.6 , f3 = 0.13 . (b) 
Relaxation to the bulk state u — u~ ; E = 0.2, (5 = 0.17. 
For both cases the system length is L = 30 , the number of 
grid-points is 300 , time step dt — 2 x 10~ J . 

The bifurcation diagram of one- and two-site nuclei 
along with a single front in the (E , (3) plane is shown in 
Fig. |^. The single kink is pinned below its bifurcation 
line and turns into a propagating burst wave above the 
line. Each nucleus is stable below its bifurcation line, 
which consists of two branches. Crossing the right or 
left branch upward results in the nucleation of a pair of 
burst waves or relaxation to the bulk state, u = w_ , 
respectively. The two branches intersect at a singular 
point, known as a cusp catastrophe 0. We see in Fig. ^ 
that the nucleation of burst waves from localized nuclei 
demonstrates hysteresis, i.e. there is a parameter range, 
where a localized nucleus is stable, while a single front 
already propagates. The bifurcation lines for multi-site 
nuclei look similar. The field value at the singular point 
on each of the lines goes to oo and (3 — > 0, with the 
number of sites in the nucleus, due to the increasing sep- 
aration of kinks. Note that the hysteresis and the nu- 
cleation phenomena mentioned above happen below the 
threshold driving strength (\E\ = 1 in Fig. |J), at which 
the effective potential from Eq. (^) loses its minima. 

For any stationary solution, the dynamic equation (||) 
reduces to the following tridiagonal system of algebraic 



(3 (u i+ i + Ui-i - 2ui) - sinuj + E 







(4) 



This arises because, in the stationary regime, Eq. (|^) 
between the sites turns into a one-dimensional Laplace 
equation with linear solutions. 

We notice that the maximum value of (3 on the one- 
site nucleus bifurcation line in Figure ||, is ~ 0.25, i.e. 
deep in the discrete (small j3) regime, which ends near 
j3 ~ 1 . It is, therefore, natural to study the stability of 
the nucleus, taking [3 as a small parameter of the analysis. 
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FIG. 2. Numerically obtained bifurcation diagram of kink 
and one- and two-site nuclei. Parameters are as in Fig. |l| 

The following three sites are the key elements of our 
analysis: the site of the nucleus, i — n , and its two near- 
est neighbors, i = nil , with the corresponding field val- 
ues u n and u n ±\ . The symmetry of the nucleus implies 
that « n+ i = u„_i = u* . The neighboring sites are both 
equivalent to the "front site" on a single kink. They are 
most active in a sense that they are the first candidates 
to cross the threshold value uq, if the stationary solu- 
tion becomes unstable. It can be shown (cf. results of 
Ref. |^] for a single kink) that the field value at the other 
sites decays with the distance from the front site down 
to u_ , as u n ±k — u- ~ (3 k ~ 1 (for k = 2, 3, . . .). Consider 
the stability analysis to first order in (3 . Then only the 
abovemcntioned three sites make a nontrivial contribu- 
tion to the system stability properties, which leaves only 
two independent variables, u n and u* . Using the corre- 
sponding two equations from the system (Q), we reduce 
them to one equation for the "active" value u* 



g{u*) = 2(3 {u- -u*) - 2sinu* + 3E 

sin 2u* — u_ H — (sin it* — E) 

P 



= 



(5) 



The solution of Eq. (|s|) is represented graphically in 
Fig. |^. The case of the right (breaking) branch of the 
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one-site nucleus bifurcation line (see Fig. ||) corresponds 
to Fig. |^(a). Then below the bifurcation line, the function 
g{u) in Eq. (||) is given by the solid line, and the equation 
has four solutions (open circles in the figure). To deter- 
mine the stability of these solutions, one has to linearize 
system (^) and find the eigenvalues of the obtained linear 
problem. We note that the maximum eigenvalue X max is 
responsible for the stability of the most active front site 
with u — u* . This means that near the bifurcation line 
only Xmax can change sign, whereas all other eigenval- 
ues remain negative. To find X m ax 7 we linearize Eq. (|5|) 
and conclude that the 1st and the 3rd solutions for u* 
have Xmax < , while the 2nd and the 4th solutions have 
Xmax > . The two former solutions are, therefore, the 
stable nodes, and the two latter solutions are saddles with 
one unstable direction in the phase space. 
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FIG. 3. Solution of Eq. (§). (a) E = 0.6 , (3 C = 0.121 ; (b) 
E = 0.2, (3 C = 0.159. 

The 1st solution, however, is trivial in the sense that 
all sites are in the bulk state u = U- . This leaves the 
single nontrivial stable solution, u nuc i in Fig. ^|. As (3 in- 
creases, the stable solution and the right saddle, u = Ub r , 
approach each other and eventually merge at (3 = J3 C . 
This occurs when the maximum of curve g(u) in Fig. |3j(a) 
touches the line g = , at u = u cr a ■ Apparently the two 
corresponding maximum eigenvalues both become zero, 
which implies a saddle-node bifurcation. Further increase 
of (3 leaves the system without a stationary solution, and 
the nucleus breaks, leading to the formation of a pair of 
oppositely moving burst waves [Fig. |l|(a)]. 

Stability analysis of the left (relaxation) branch of the 



bifurcation line (see Fig. |J) can be made analogously, 
using Fig. ||(b). The difference is that now the stable 
solution merges with the left (relaxation) saddle, u = 
u re i , and the nucleus relaxes to the bulk state [Fig. 0(b)]. 

It appears that the first-order analysis is insufficient for 
the quantitative comparison of our theory with numerical 
simulations of Eq. (Q). To improve our predictions, we 
have developed the stability analysis to second order in 
[3 . In this case one has to consider the nucleus and two of 
its neighbors from each side. Symmetry of the problem 
reduces the number of independent variables from five 
to three, and we deal with three corresponding equations 
from system (||) . We then perform the same procedure as 
for the first-order stability analysis described above. In 
Fig. [| we plot the resulting bifurcation line for the single- 
site nucleus along with the results of the simulations from 
Fig. @. We see in Fig. || that our theoretical prediction is 
in good quantitative agreement with the simulations. 
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FIG. 4. Comparison between analytical and numerical bi- 
furcation lines for a one-site nucleus. Parameters as in Fig. 

We now describe the physical mechanisms of the ob- 
tained phenomena of relaxation and break-up of local- 
ized nuclei. In the parameter regime below the bifur- 
cation line, the energy functional £, Eq. ([!]), possesses 
a minimum in the functional space, which corresponds 
to a stable nucleus. Separatrices connect this minimum 
to the two neighboring unstable solutions (saddles). If 
one starts with an initial condition in cither of these 
saddles slightly perturbed in the direction of the min- 
imum, the system evolves along the separatrix until it 
relaxes at the minimum. If one parameterizes the po- 
sition on the separatrix with an arclength s, defined as 
ds 2 = dt 2 J(du/dt) 2 dx , the system dynamics take a sim- 
ple gradient form, ds/dt = —d£(s)/ds . 

In Fig. U we plot the dependence of energy E on s 
from numerical simulations. Fig. [5](a) shows £ (s) near 
the break-up branch of the bifurcation line. We see in 
the figure that, for the stable nucleus (solid line), the en- 
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ergy indeed has a minimum and two neighboring maxima, 
corresponding to the saddles of the dynamics. When the 
nucleus loses stability (dashed line) , the minimum merges 
with the right (break-up) saddle. Then the nucleus turns 
into a pair of burst waves and the energy begins mono- 
tonically decreasing. If we start with an initial condition 
slightly to the right of the break-up saddle then the sys- 
tem develops a pair of burst waves, even in the stable nu- 
cleus regime. This is the manifestation of the nucleation 
hysteresis effect, since the bound nucleus is still stable 
in the parameter range for which a well-separated front 
already bursts and propagates. If, on the other hand, we 
start to the left of the relaxation saddle, then the system 
will relax to the homogeneous bulk state u = it_ for both 
stable and unstable nuclei. 



0113<P C 
0.13>p 



relaxation 



CO -5 



-10 




5 10 
S 



15 



p=o.i<p c 

P=0.2>p c 



relaxation 



3— site nucleus 



2 - 




10 



15 



FIG. 5. Energy £ vs arclength s near (a) break-up branch 
(E = 0.6) and (b) relaxation branch (E = 0.2) of the bifur- 
cation line of a one-site nucleus. 

The system behavior near the relaxation branch of the 
bifurcation line is different, see Fig. ||(b). As crosses 
the critical value, the minimum of the functional merges 
with the left (relaxation) saddle, and the system relaxes 
to the homogeneous bulk state. The same thing happens 
if we start to the left of this saddle, for the regime of 
parameters where the nucleus is stable. However, if we 



start slightly to the right of the break-up saddle, then the 
nucleus at first starts breaking, but then relaxes to the 
three-site nucleus, rather than developing into the pair 
of separated burst waves. 

In conclusion, we have demonstrated that nonlinear 
lattices with diffusive coupling possess localized struc- 
tures, in certain parameter regimes. We have shown that 
these structures can be destroyed in two alternative sce- 
naria: (i) break-up into a pair of oppositely propagating 
burst waves, and (ii) relaxation to a homogeneous state. 
We have found hysteresis in the burst-wave nucleation 
from a localized embryo, appearing as the difference be- 
tween the stability thresholds of nucleus break-up and 
burst-wave propagation. We have given a theoretical de- 
scription of these phenomena in terms of the energy func- 
tional of our system. The predictions of our theory are 
in a good quantitative agreement with numerical simula- 
tions of the full system. Our interesting directions for fu- 
ture research: the influence of the localized structures on 
thermal nucleation [||, dynamics of localized structures 
in corresponding two- and three-dimensional systems, ex- 
perimental studies of stable and unstable nuclei in ma- 
terials, etc. We note that, though we have studied lat- 
tices with continuous coupling, the stability analysis has 
been performed for stationary localized structures from 
Eq. (|J). Therefore, the obtained stability properties of 
these structures should remain intact for completely dis- 
crete systems, such as arrays of Josephson junctions Q 
or lattices of nonlinear oscillators Q . 
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